Thermodynamic properties of the itinerant-boson ferromagnet 
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Thermodynamics of a spin-1 Bose gas with ferromagnetic interactions are investigated via the 
mean-field theory. It is apparently shown in the specific heat curve that the system undergoes two 
phase transitions, the ferromagnetic transition and the Bose-Einstein condensation, with the Curie 
point above the condensation temperature. Above the Curie point, the susceptibility fits the Curie- 
Weiss law perfectly. At a fixed temperature, the reciprocal susceptibility is also in a good linear 
relationship with the ferromagnetic interaction. 

PACS numbers: 75.40.Cx, 75.10.Lp, 75.30.Kz, 03.75.Mn 



I. INTRODUCTION 

The realization of spinor Bose-Einstein conden- 
sation in optical traps^ 2 - has stimulated enormous 
interest in magnetic properties of quantum Bose 
gaS es3^&LS^iLi2^^. In optical traps, the hy- 
perfine degree of freedom of confined atoms, such as 
87 Rb, is released and therefore the atom can exhibit mag- 
netism. More intriguingly, an exchange-like spin-spin in- 
teraction can be present between atoms. In the F = 1 
87 Rb atoms, the interaction is ferromagnetic^, so the 
87 Rb gas appears to be a prototype of itinerant-boson 
ferromagnet^^i 7 ^. 

Ferromagnetism is one of the central research themes 
in condensed matter physic a 16 ' 17 . Two types of ferro- 
magnetism have already been intensively studied: local- 
moment ferromagnetism and itinerant-electron ferromag- 
netism. Although particles in these two systems obey 
different statistics, they both share some common fea- 
tures. For example, both ferromagnets have a Curie 
point, above which the susceptibility conforms to Curie- 
Weiss law. Nonetheless, from the theoretical point of 
view, the origin of Curie- Weiss law is quite different 
for these two systems. In insulators it is due to local 
thermal spin fluctuations and can be easily explained in 
the mean-field approximation. On the other hand, in 
itinerant-electron ferromagnets the Curie- Weiss law may 
be caused by the mode-mode coupling between spin fluc- 
tuations and the theoretical treatment is much more com- 
plicated 17 . An appropriate theory is the self-consistent 
renormalization (SCR) theory 18 which goes beyond the 
Hartree-Fock approximation and the random-phase ap- 
proximation. The SCR theory succeeds in explaining var- 
ious magnetic properties of itinerant-electron ferromag- 
nets and is also extended to treat the specific heal^. 

The 87 Rb gas provides opportunity to study the third 
type of ferromagnetism. Ho^, Ohmi and Machida 5 have 
studied its ground state properties and the spin-wave 
spectrum. The long wavelength spectrum is linear in k, 
the wave vector, as in the two former cases. In our previ- 
ous papers, we have investigated the finite-temperature 
properties, especially the Curie poim— . We suggest that 
the phase diagram in itinerant bosons should be more 
complicated than the other two ferromagnets, because 



the Bose system has an intrinsic phase transition, other 
than the ferromagnetic transition. An interesting conclu- 
sion we arrived is that its Curie point, Tp, is never be- 
low the Bose-Einstein condensation temperature, Tc, re- 
gardless of the magnitude of the ferromagnetic coupling 8 . 
Kis-Szabo el al got the same point later-. However, ther- 
modynamics of the itinerant-boson ferromagnet has not 
yet been investigated systematically so far. 

The purpose of this paper is to calculate the thermo- 
dynamic quantities of ferromagnetic bosons. As in the 
fermion case, the specific heat and magnetic susceptibil- 
ity are of the most interest. In Section 2, we introduce 
the mean-field approximation to deal with ferromagnetic 
interaction, taking the spin-1 Bose gas as an example. 
In Section 3, phase transitions are discussed by calcu- 
lating the free energy and specific heat. In Section 4, 
the susceptibility above the Curie point is calculated. A 
summary is given in the last section. 



II. THE MEAN-FIELD APPROXIMATION 

The spin-1 Bose gas with ferromagnetic couplings is 
described by the following Hamiltonian, 



h e Vv(r) 



-~I S I drS(r) • S(r), 



(1) 



where i/v(r) is the quantum field operator for annihi- 
lating an atom in spin state \a) at site r. For a spin-1 
gas, a = +1,0,-1. The parameter h e denotes the ex- 
ternal magnetic field. The last term represents the ferro- 
magnetic exchange between two different bosons meet- 
ing at site r and I s (> 0) is the exchange constant. 
S = {S x , S v , S z } are the spin operators, which can be 
expressed via the 3x3 Pauli matrices, for example, 
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Within the mean-field approximation, we treat the spin- 
dependent interactions as a molecular field except of a 
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particle with itself, 

- -S- S Ps -IS) • s 



(S) = -MS Z + ^M 2 , (3) 



where M = (S z ) is the ferromagnetic order parameter. 
Then the effective Hamiltonian for the grand canonical 
ensemble reads, 

H - Nu = ^[e k - u - a(h m + h e )]n ka + ^M 2 I S N, (4) 

where ek is the kinetic energy for free particles, h m — 
I S M is called the molecular field, similar to the Stoner 
theory for fermion gasesi£; n is the chemical potential; N 
is the total particle number. The grand thermodynamic 
potential can be worked out in a standard way, 
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where h = h rn + h e , m is the mass of particle, and / is 
the polylogarithm function defined by 



fn(x) = J2 



x \ k 
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(6) 



where x < 0. The mean- field self-consistent equations 
are derived from the grand thermodynamic potential, 
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(7a) 



(7b) 



where n is the density of particles, no is the density of 
condensed one and M = ±iM- is the magnetization, no is 
zero unless the temperature is below the BEC point T c . 



III. THE FREE ENERGY AND SPECIFIC HEAT 

In our previous investigations, we showed that the sys- 
tem exhibits two phase transitions, the Bose-Einstein 
condensation (BEC) and the ferromagnetic transition 8 . 
The condensation temperature Tc and the Curie tem- 
perature Tp are calculated by solving the self-consistent 




FIG. 1: Free energies of the FM and PM states with the 
ferromagnetic coupling 1 = 1. The two curves cross at the 
temperature tf ft! 0.80, which is just the FM transition point. 



equations. We find that Tp is never below Tc for all 
systems with a finite ferromagnetic exchange {I s ^ 0). 

However, one can get another solution to the Eqs. (J7J), 
with M = at all temperatures. It means the system 
does not undergo a ferromagnetic transition at all, but 
remains in paramagnetic (PM) state at low temperatures. 
Actually, whether there exists a Curie point in the fer- 
romagnetic Bose gases is still a controversial question. 
Some researchers suppose that the Bose gas can not be 
magnetized spontaneously at low temperatures even if 
the ferromagnetic exchange is present^. 

In order to single out the physically correct solution, 
one has to compare the free energy of the the ferromag- 
netic (FM) state and the PM state. The relation between 
the free energy and the grand thermodynamic potential 
has the form: 



f = n + Nu 



(8) 



For computational convenience, the temperature T and 
exchange interaction I s are re-scaled, as did in Ref. [8], 

by the following formula: t = [3C(|)] _ ^T/T and I = 
[3C(|)r f VfeTo), where 
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is the condensation temperature of ideal spin-1 Bose gas. 
Hereinafter, all the numerical results are obtained by set- 
ting n = k B = = 1. Figure 1 shows the free energy 
of unit volume for the gas with / = 1.0. It shows clearly 
that the free energy of FM state is lower than that of the 
PM state at the low temperature region, which demon- 
strates that the FM state should be more stable than 
PM state. Therefore, the low temperature state has a 
spontaneous magnetization. In experiments, the total 
spin of the ferromagnetic spinor condensate is observed 
to be conserved, which is called the spin conservation rule 
in some literature o 13 i 14 . However, the spin conservation 
rule holds only globally, not locally. In the theoretical 
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FIG. 2: Specific heats of spinor Bose gases with the coupiing 
J = 1.0 and 0.4. The dotted vertical lines serve to guide the 
eye to see the transition points. 



treatment of Ref . [20] , the spin conservation rule is im- 
posed by introducing a lagrangian multiplier. It is over- 
constrained in some sense, so that the spontaneous mag- 
netization can not be established. Recent experiments 
and theories indicate some domain structures should be 
formed and each domain is magnetize d 11 ' 12 ! 15 , where the 
conservation law for the total spin can be restored natu- 
rally. 

The FM transition is induced by the FM coupling and 
the transition temperature is about tf ~ 0.8 for the Bose 
gas with / = 1.0. When the temperature goes down fur- 
ther, the BEC then occurs, which is the intrinsic phase 
transition of Bose gases. To demonstrate different fea- 
tures of the two transitions, we now calculate the specific 
heat of unit volume, 
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where U is the internal energy 
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As shown in Fig. 2, for the system with / = 1.0, the 
specific heat exhibits a jump discontinuity at tf w 0.8, 
from the PM state to the FM state. This is a charac- 
teristic feature of the Landau-type of second-order phase 
transition. And similar behaviors have been observed in 
the specific heat of ferromagnetic insulators or itinerant - 
fermion ferromagnet a 16 ' 17 . The BEC occurs at t c » 0.5, 
where the specific exhibits a bend. But specific heat is 
continuous at the BEC point, similar to that of a free 
Bose gas. The results indicate that the critical behaviors 
are different at the two transition points on the mean- 
field level. 



IV. THE CURIE- WEISS LAW 

For a ferromagnet, the susceptibility above the Curie 
point is of special interest. As already studied, the sus- 
ceptibility is well described by Curie- Weiss law both in 
the insulating ferromagnet and the itinerant-electron fer- 
romagnet. In this section we calculate the susceptibility 
for the itinerant-boson ferromagnet. 

The susceptibility can be derived from Eqs. ([7]). Dif- 
ferentiating both sides of the two equations and removing 
the term dfi, the following equation are deduced, 



dM = fd 



I S M + h e 
k B T 



(11) 



where 

/ = 
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Above the Curie point, the magnetization M (then h — 
I S M + h e ) diminishes correspondingly when the external 
field h e tends to zero. So the second term in the above 
equation is omitted and then / has a simple form: 
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(13) 



Thus the zero-field susceptibility of unit volume is given 
by 



( dM 
\dh~e 



1 
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(14) 



The susceptibility % is a function of the coupling I s and 
temperature T. Figure 3 shows 1/x and \ versus I at dif- 
ferent given temperatures. As shown in the inset of Fig. 

3, the susceptibility becomes larger as the coupling I in- 
creasing. It is physically reasonable since the the system 
with larger / can be magnetized more easily. At a given 
temperature, \ diverges as / approaches a critical value. 
It is worth noting that the inverse of the susceptibility is 
in a good linear relationship with the coupling. 

The susceptibility versus temperature is shown in Fig. 

4. One can immediately find that the susceptibility meet 
quite well with the Curie- Weiss law in a very large tem- 
perature region. Seeing that the Curie- Weiss law is very 
difficult to be derived for the itinerant-fermion ferromag- 
net, it is really surprising that we get it for the itinerant- 
boson ferromagnet just based on the mean-field approx- 
imation. 

In order to discuss the Curie- Weiss law in a more ex- 
plicit way, we proceed to carry out a semi-analytical cal- 
culation to deduce the linear dependence of 1/x on the 
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FIG. 3: Magnetic susceptibilities versus ferromagnetic cou- 
plings of spinor Bose gases at temperature t = 0.4, 0.6 and 
0.8. 




FIG. 4: Magnetic susceptibilities versus temperatures of 
spinor Bose gases on coupling I = 0.2, 0.4 and 1.0. 



temperature. The first step is to analyze the tempera- 
ture dependence of /. It is quite complicated, because the 
chemical potential /i is an implicit function of the tem- 
perature. We consider a limit case that the parameter I s 
is quite small, when Tf is close to T c . So [i is close to zero 
in the vicinity of Tf . According to the asymptotic behav- 
ior of the polylogarithm function: fs(x) « C(|) — 2^/nx 

and fi(x) rs y/n/x as x — » CP, we get the following 
equations from Eqs. (7a) and (13) respectively, 
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and 



t ' k B Tm \ 
- ' 2nh 2 J 



ksTir 



(16) 



Substitute Eqs. HU) and into Eq. HI]), we get 
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In the vicinity of Tp which is only slightly larger than 
Tq, Eq. (fTTj) could be further simplified to 
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Thus the effective FM transition temperature is defined 
as 



T f =T n 



8tt 



9c 2 (|)fcB 



Is- 



So far the Curie- Weiss law is derived. We note that the 
derivation is only valid in small I s cases. 

In the high temperature limit, one can also easily prove 
that x _1 is linearly dependant on T. In this case, — r^j? 
has a quite large value, so that 



A 4 



5 V k R T 



li 



5 V k R T 



according to Eq. (6). Combining Eqs. (7a), (13) and 
(14), it yields 

X - l =rr 1 {k B T-I s ). (19) 
We estimate this equation holds in the range of t > 10. 

V. SUMMARY 

In summary, we calculate thermodynamic quantities 
of the spinor Bose gas with ferromagnetic interactions. 
Such kind of investigations has already been performed 
intensively for the ferromagnetic fermions, while few as 
yet for bosons. Based on a mean-field approximation, 
we show that the system undergos a ferromagnetic phase 
transition first, then the Bose-Einstein condensation with 
the temperature decreasing. The specific heat shows a 
jump discontinuity at the Curie point and a bend at the 
Bose-Einstein condensation temperature, indicating that 
critical behaviors are different near the two transition. 
The more surprising result is that the mean-field theory 
yield the magnetic susceptibility which satisfies perfectly 
the Curie- Weiss law over a wide range of temperature. 
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